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Abstract. Observations of the X-ray pulsar 4U 2206+54, performed over a 15-year span, show its 
period, which is now 5555 ± 9 s, to increase drastically. Such behavior of the pulsar can be hardly 
explained in the frame of traditional scenarios of spin evolution of compact stars. The observed 
spin-down rate of the neutron star in 4U 2206+54 proves to be in good agreement with the value 
expected within magnetic accretion scenario which takes into account that magnetic field of the 
accretion stream can under certain conditions affect its geometry and type of flow. The neutron star 
Q_j in this case accretes material from a dense gaseous slab with small angular momentum which is 

kept in equilibrium by the magnetic field of the flow itself. The magnetic accretion scenario can be 
realized in 4U 2206+54 provided the magnetic field strength on the surface of the optical counterpart 
to the neutron star is in excess of 70 G. The magnetic field strength on the surface of the neutron 
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star in the frame of this scenario turns out to be 4 x 10 G which is in accordance with estimates of 



this parameter from the analysis of X-ray spectra of the pulsar. 
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g: 1. INTRODUCTION 

4U 2206+54 is a moderate-luminosity X-ray source discovered at the beginning of X- 
ray astronomy [1]. It has been identified with a neutron star in a high-mass binary with 
the orbital period P or b ~ 19.25 d [3] situated at the distance of about 2.6 kpc [2]. Its 

o 

optical counterpart is the 09.5V star BD +53 2790 underfilling its Roche lobe and 
^ ; losing mass via relatively slow, v w ~ 300 — 500kms _1 , dense stellar wind [0, \$\. The 

neutron star rotates with the period P s = 5555 ± 9 s [@, 0]. The X-ray spectrum of the 
system is typical for a neutron star accreting material onto its magnetic poles. The 
average X-ray luminosity throughout the entire duration of observations ranges within 
L x ~ (2 — 4) x 10 35 ergs _1 . The individual X-ray light curves exhibit chaotic variations 
of relatively low amplitude on a timescale of 30 minutes. Higher-amplitude fluctuations 
of 15 to 45 hours duration and a peak luminosity of ~ 4 x 10 36 ergs~ ! were detected in 
hard X-rays with INTEGRAL mission H]. 

Analyzing X-ray light-curves of 4U 2206+54 obtained in 1998-2007 Finger et al. [0] 
reported the spin frequency of the neutron star (v = 1 /P s ) to decrease at the rate v s d ~ 
( — 1.7±0.7) x 10 _14 Hzs _1 . The characteristic spin-down time-scale of the neutron star, 
T S( j = P S /2P ~ 180 yr, in this case is very short compared to a typical age of neutron stars 
in the high-mass X-ray binaries (> 10 6 - 10 7 yr, see, e.g. |9(]). This testifies the episodic 



character of observed variations of the spin period and allows to assume that the star 
rotates near its equilibrium period. It is also worth mentioning that the X-ray luminosity 
of the system during this time remained close to its average value with except for hard X- 
ray flares detected at the end of 2005 [8]. This indicates an absence of direct correlation 
between variations of the star's spin period and mass accretion rate onto its surface 
and marks it possible to discard a possibility to describe the observed spin evolution 
exclusively in terms of variations of the spin-up torque applied to the star from the 
accretion flow (for discussion see [0]). 

Interpretation of the observed spin-down rate of the neutron star in the frame of 
popular accretion scenarios was discussed in [0, Oil]. In particular, the authors argued 
that a presence of the Keplerian accretion disk is the system is hardly probable. The 
spin-down of the star in this case could occur only under condition that the magnetic 
field strength on its surface exceeds 10 15 G. In order to explain the observed spin- 
down rate of the star within the picture of free-falling quasi- spherical accretion flow 
it is necessary to assume that the value of the surface magnetic field of the neutron star 
is in excess of 10 14 G. Finally, it was shown that even within the model of accretion 
from a turbulent spherical envelope the observed behavior of the pulsar can be described 
only provided the magnetic field strength on the surface of the compact star is at least 
6 x 10 13 G and, thus, surpasses the quantum critical limit to the magnetic field strength 
Bq = mlc 3 /eh ~ 4.4 x 10 13 G. The neutron star in the frame of this scenario turns out to 
be an accreting magnetar of the age ^ 30000yr estimated from the characteristic decay 
time of supercritical magnetic field 

An assumption about supercritical value of magnetic field in neutron stars seems quite 
plausible in the light of investigation of anomalous X-ray pulsars and soft gamma-ray 
repeaters. However, it is not obvious that this assumption is acceptable in the case of 
long -period X-ray pulsars. In particular, analysis of the X-ray spectrum of 4U 2206+54, 
presented in 111 ill , has shown that surface magnetic field of the neutron star in this 
system is Bqrsf — 4 x 10 12 G. This is testified by the observed energy distribution in 
the spectrum of this source and a presence of absorbtion feature in the region of 30keV 
identified by the authors with a cyclotron line. It should be noted that in case of relatively 
slow stellar wind, v w ~ 300 — 500km s _1 , and low rate of mass exchange between the 
system components, 

M a ~2xlO l5 mR 6 ( if T ) gs- 1 , (1) 
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the spin-down time-scale of the star in the propeller state (see equation (21) in 1120), 



exceeds the decay timescale of supercritical magnetic field (see |fl3l IbH]). Here m, R$ 
and V45 are the mass, radius and the moment of inertia of the neutron star expressed in 
units of 1.4Mq, 10 6 cm and 10 45 gcm 2 , respectively, while /I32 is the dipole magnetic 
moment of the star in units of 10 32 Gem 3 . M15 = 10 15 g s _1 is the rate of mass exchange 
between the system components which in the case of stationary accretion equals the 



mass accretion rate onto the neutron star surface. This result raises doubts concerning 
the validity of the assumption about extremely high magnetic field of this object. 

A situation in which the magnetic field strength obtained in the spin evolution models 
significantly exceeds the value estimated from the cyclotron line measurements is not 
unique. For instance, such discrepancy was revealed for a long-period X-ray pulsar 
GX 301-2 fli. Our studies [180 have shown that this inconsistency can be 
caused by oversimplification of the accretion picture used in the modeling of spin 
evolution of X-ray pulsars in massive binary systems. In particular, this flaw can be 
avoided in the frame of magnetic accretion scenario [fl9l Eol 21] accounting for the 
influence of the intrinsic magnetic field of the accretion flow on its geometry and 
parameters. In this paper we show that conditions necessary for the magnetic accretion to 
be realized in 4U 2206+54 can be satisfied if the magnetic field strength on the surface of 
the massive component exceeds 70 G. The spin-down rate of the neutron star evaluated 
within this scenario corresponds to its observed value in case the magnetic field strength 
of the neutron star is of the order of ~ #crsf- 



2. MAGNETIC ACCRETION IN 4U 2206+54 

4U 2206+54 is a close binary system with wind-fed accretion. The neutron star which 
moves through the stellar wind of its massive companion with the relative velocity v re i, 
captures matter at the rate M = nE^p^y^. Here Rq = 2GM ns /v^ el is the gravitational 
(Bondi) radius of the neutron star, M ns is its mass and is the wind density at the 
Bondi radius. As the captured material is approaching the neutron star under the force 
of its gravitational attraction it interacts with the stellar magnetic field. This leads to 
the formation of magnetosphere of radius r m which is defined by equating the ram 
pressure of the accretion flow with the magnetic pressure in the neutron star's dipole 
field. The magnetosphere in the first approximation prevents the accretion flow from 
reaching the stellar surface. The flow is decelerating at the magneto spheric boundary 
and after penetrating the magneto spheric field and moving along the field lines reaches 
the neutron star surface in the region of magnetic poles. 

The radius of magnetosphere depends on the structure of accretion flow beyond its 
boundary which is determined by the following key parameters: the sound speed, c s q, 
and the magnetic field strength BfQ, in the matter captured by the neutron star at the 
Bondi radius, as well as the relative velocity of the neutron star with respect to the 
stellar wind of its massive counterpart, v re \, and the orbital angular velocity of the binary 
system Q or b = 2n/P OI \). Depending on relationship between these parameters the process 
of accretion proceeds either via a Keplerian disk formation or in the quasi- spherical or 
magnetic regime. 



2.1. Non-magnetic accretion 

In the conventional (non-magnetic) accretion scenarios the magnetic field in the 
matter captured by the neutron star at the Bondi radius is assumed to be negligible 
and not to influence the structure of accretion flow. The model of spherical accretion 



within this approach (a so called Bondi accretion) was first introduced in 11221 12311 and 



was further used in the modeling of accretion in the X-ray pulsars [24, 25]. In the frame 
of this model the accreting material is approaching a neutron star in a state of free- 
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fall. Moving at the free-fall velocity Vff = (2GM ns /r) 1 it reaches the surface of the 

1 /2 

neutron star at dynamical (free-fall) time-scale tff = (r 3 /2M ns ) . The ram pressure in 
the accreting flow is increasing as it approaches the neutron star as [12311 

[R \ 5/2 

E tam {r)=E tam {R G )l-±) , (3) 

where E ram (Ro) = pooV 2 el is the ram pressure at the Bondi radius. The balance between 
the ram pressure of the spherical flow and the magnetic pressure of the neutron star 
dipole field P m = ^/Inr 6 , is reached at the radius [12411 
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where /i is the dipole magnetic moment of the neutron star. This is a magneto spheric 
radius of the neutron star in the non-magnetic accretion approximation. 

Quasi- spherical accretion is a modification of Bondi accretion accounting for a non- 
zero angular momentum of matter captured by a neutron star in a binary system. The 
process of mass accretion in this case is accompanied by the angular momentum accre- 
tion at the rate J = t, Qo^RqM (see ll26ll and references therein), where ^ is a parameter 
accounting for angular momentum dissipation in the accretion flow. The average value 
of this parameter calculated in the numerical modeling of the wind-fed accretion un- 



der assumption about a zero magnetic field of the flow is < <£; >= 0.2 (see Q2VD and 
references therein). 

Deviation of the quasi-spherical flow geometry from spherical symmetry remains 
insignificant under condition r c i rc < r a , where r c i rc = J 2 /GM ns M 2 is a circularization 
radius, at which the angular velocity of the accreting material, co en = <r ^orb (Rc/ r ) > 

1/2 

equals the Keplerian angular velocity, ft^ = (r 3 /2GM ns ) . Otherwise, a Keplerian 
accretion disk forms in the system. The condition of disk formation (r c j rc > r a ) can be 
written as v re i < v cr , where Ill2n 



v cr „ 210 „ 3 -V» „!./» Mf » Qt) ^ kms -.. (5) 

Here (^0.2 = <d/0-2. For the parameters of 4U 2206+54 we get the value of v cr signif- 
icantly less than the stellar wind velocity derived from observations. This means that 
formation of the Keplerian accretion disk is rather unlikely from the theoretical point of 
view. In this light an absence of observational manifestations of the Keplerian disk in 
the system seems rather natural flUl . 



2.2. Magnetic accretion 



As first shown by Shvartsman 111911 . the picture of wind-fed accretion in the high-mass 
X-ray binaries can essentially differ from expected within Bondi scenario if the accreting 
material is magnetized. In the free-falling accretion flow magnetic field is predominantly 
radial. This is caused by the fact that in the process of spherical accretion all transverse 
scales contract as ~ r~ 2 while the longitudinal scales expand as~ r 1 / 2 ll28ll . Then under 
the condition of mag netic flux conservation the field strength in the accreting matter 
rises as B r r~ 2 [29]. This means that the magnetic energy density E m = B^/Sn, in the 
free-falling plasma, 

E m (r)=E m {R G )(—) , (6) 
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increases more rapidly than the kinetic energy of its radial motion E mm {r) °c (R G /r) 1 
(see expression©. Here E m (Ro) = j5^ l E t h(RQ) is the magnetic energy density in the 
accretion flow at the Bondi radius normalized using the parameter /3 to the thermal 
energy of the flow, E±(Rq) = Pc»c 2 . The distance at which the magnetic energy density 
in the accretion flow balances its kinetic energy density (the Shvartsman radius) can be 
obtained solving the equation E m (R s h) = E iam (R s h), in the form 

_ fl -2/3 2GM m c% 3 

R sh~P 1073—- ( 7 ) 
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Shvartsman | fl9]| has noted that mass accretion in the spatial region r < R s h can take 
place as long as the magnetic field in the flow dissipates. Otherwise the magnetic energy 
in the accretion flow would exceed its total energy (equal to gravitational energy of the 
star at a given radius) which is impossible under the energy conservation law. Thus, the 
velocity of radial motion in this region cannot exceed v r ~ r/t KC , where 

77m 1 hi (^] (8) 



is the time of field dissipation due to magnetic reconnection, and va = B t / (ATZp) 1 / 2 is 
the Alfven velocity in the accretion flow. T] m is a parameter allowing for the efficiency 
of reconnection. Its value depends on both the physical parameters in plasma and 
configuration of the magnetic field and in general case ranges as < T} m < 0.1 ||30)l . 
Since va < Vff (the equality is met at the Shvartsman radius) the time of field dissipation 
remains significantly less than dynamical timescale, tg, throughout the accretion process. 
On one hand this justifies an assumption about the magnetic flux conservation in the 
free-falling plasma, but on the other hand this means that accretion flow is decelerated 
by its own magnetic field at the Shvartsman radius and further matter inflow proceeds in 
the diffusion regime. 

Basic conclusions of the scenario proposed by Shvartsman were later confirmed in 
quantitative modeling by Bisnovatyi-Kogan and Ruzmaikin [2(J 2j]], and by the results 



of numerical simulations of the spherical magnetic accretion presented in ||3 lL 13211. These 



authors have shown that magnetic field amplification in the free-falling material leads 
to deceleration of the accretion flow and its shock-heating up to adiabatic temperature. 
Accretion inside the Shvartsman radius occurs on the timescale of magnetic field dissi- 
pation. The accretion picture in this case depends on the efficiency of cooling processes 
in the accreting material. If the cooling time at the Shvartsman radius, ? C ooi(^sh)> ex- 
ceeds the heating time due to magnetic energy dissipation, t KC (R s h), the accretion flow 
transforms into a hot turbulent envelope with some portion of matter leaving the system 
in a form of jets |3ll [32 ]. Otherwise, the flow forms the magnetic slab whose struc- 



ture is determined by geometry of the large-scale magnetic field in the accreting matter 



112a 1210. 

In the X-ray pulsars the magnetic field of the accretion flow influences its structure 
in case > r a . This condition is met if the relative velocity of the neutron star with 



respect to the wind of its massive companion satisfies inequality v re i < v mca , where [|18l1 

v mca = r 1 / 5 (2GM ns ) 12 / 3 V 6 / 35 M 3 / 35 c 2 ' 5 . (9) 
Substituting parameters of 4U 2206+54, we find 

v™ * 460 r 1 ' 5 <"' 2/3 W 35 M?f doi^) 2 ' 5 km.-'. (10) 

Thus, the accretion process in this system would occur according to magnetic accretion 
scenario if /3 ~ 1 . Estimating the thermal energy density of the matter captured by the 
neutron star at the Bondi radius, 

and its magnetic energy density, 

*.(«) = 10-3 ergo™- (^-,) 2 (^-f, (12) 

we can conclude that the accretion process in 4U 2206+54 should be described in terms 
of magnetic accretion if the strength of the dipole magnetic field on the surface of the 
massive (optical) component of this system exceeds 70 G. Here a\j, = a/10 13 cm is the 
orbital separation, /i ms is the dipole magnetic moment of the massive star and a\ is the 
distance (measured from the massive component) at which the kinetic energy density of 
the stellar wind reaches the magnetic energy density of its dipole field (for discussion see 
rfl8ll '). The free-fall of matter captured by the neutron star in this case will be terminated 
by the intrinsic magnetic field of the accretion flow at the Shvartsman radius exceeding 
the magneto spheric radius of the star. The accretion process inside the Shvartsman radius 
will take place on the timescale of magnetic field dissipation ~ t mc . 

In order to understand the structure of accretion flow inside the Shvartsman radius it is 
necessary to take into consideration that the time of its cooling due to inverse Compton 
scattering of the X-ray photons emitted from the neutron star surface on hot electrons of 
plasma surrounding its magnetosphere Q25Q , 

3nr 2 m e c 2 

tc{r) = ^Z5T' (13) 



turns out to be less than the heating time due to magnetic field dissipation in case 
L x > L CI , where 

L C ,^ 3 x 10* „./> gffi ( J|_) ergs-'. (14) 

Here m e is the electron mass, and Or is the Thompson cross-section. This implies that 
the process of magnetic accretion in 4U 2206+54 will take place according to scenario 
worked out by Bisnovatyi-Kogan and Ruzmaikin B2Q. 12111 . if r\ m ^ 0.1(r a /i? s h) 1//2 . The 
neutron star magnetosphere in this case will be surrounded by the magnetic slab in 
which the radial motion will proceed as long as the magnetic field of the accretion flow 
dissipates. The density of plasma at the inner radius of the slab in the region of its 
interaction with the neutron star magnetosphere, p s \, can be evaluated from the balance 
of magnetic pressure and the thermal pressure in the slab as 

u 2 m D 

= oJl T 6 • < 15 ) 
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Here m p and k# are the proton mass and Boltzmann constant, and Tq is the gas tempera- 
ture at the inner radius of the magnetic slab. 



3. SPIN EVOLUTION OF THE PULSAR 

One of the signs of magnetic accretion in 4U 2206+54 is extremely high value of 
angular momentum dissipation in the accretion flow. This conclusion is made taking into 
account that spin-down of neutron stars undergoing quasi- spherical accretion can occur 
provided the angular velocity of accreting material at the magneto spheric boundary, 



O)em(fm) = £^orb C^cAm) 2 * is less than the angular velocity of the star, co s (see ll33ll ). 
Solving inequality (O em (r m ) < co s for we find 



« < O.OUn-^L^Rr x (16) 
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The derived value of £ is at least an order of magnitude less than the average value of 
this parameter evaluated in numerical simulations of the wind-fed accretion based on 
assumption /3 3> 1 (see [|27|l and references therein). So significant dissipation of the 
angular momentum in the accretion flow implies that either the neutron star is accreting 
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material from the hot turbulent quasi-static envelope [|34l [10] , or the infalling matter 
possesses hig h enough magnetic field and looses angular momentum due to magnetic 
viscosity 



The spin-up torque, applied to the neutron star from the accretion flow, ^T su 
$£l olh R 2 G M, 



K su & 10 31 ^ .0im 2 M 15 ( ^) (iro^r 4 d y nec ^ < 17 > 



under these conditions turns out to be significantly less than the absolute value of spin- 
down torque evaluated from observations, 



K sd ^ 2nlv sd ~ 6 x 10 31 / 45 ( 2 x J-u^g-i ) dynecm. (18) 

This allows to consider the equation of spin evolution of an accreting neutron star, 
16) = K su — K sd , in a simplified form 16) ~ K sd . 

The spin-down torque applied to the neutron star from the accretion flow can be 

expressed as = hM c (0^. Substituting M c = /i 2 / (2GM ns r m ) ^ 2 we get |[l8ll 

^- fr r )3/2' (19) 
V m ' cor J 
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where r cor = (GM ns /ft)^) is the corotation radius of the neutron star. The value of 
spin-down torque essentially depends on the magneto spheric radius. In case of accretion 



from a turbulent envelope the magneto spheric radius is limited as r m > r a 11241 13411 . 
and this estimate of the neutron star spin-down rate proves to be insufficient to explain 
the observed spin-down rate of the pulsar 4U 2206+54 ifioll . In case of accretion from 
the magnetic slab the value of magneto spheric radius depends on the mechanism by 
which the accretion flow penetrates into the stellar magnetic field. If the process of 
penetration is governed by interchange instabilities of the magneto spheric boundary 
(Rayleigh-Tailor and Kelvin-Helmholtz instabilities), then, as in the previous case, the 
magneto spheric radius will be close to its canonical value, r a . Otherwise it proves to 
be significantly less. If the interchange instabilities of the magneto spheric boundary are 
suppressed, the plasma penetration into the stellar field takes place due to diffusion and 
magnetic reconnection at the rate [37] 



Min(r m ) = ATZr m d m p Q Vff{r m ) = Anr 3 ^ D l £ p^' 2 {r m ) . (20) 
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Here 5m = (D e ff T d) is me thickness of diffusion layer on the magneto spheric bound- 
ary (magnetopause), D e ff is the effective diffusion coefficient and po is the plasma den- 
sity at the boundary. The characteristic time of plasma diffusion into the magnetic field, 
Td, is determined by the time necessary for plasma having penetrated into the field to 
leave the magnetopause moving along the magneto spheric field lines in the gravitational 
field of the neutron star. Under the conditions of interest this time is close to dynamic 
(free-fall) time at the magneto spheric boundary, i.e. Td ~ tff(r m ) |[3~7ll . 

The radius of the neutron star magnetosphere in this case can be evaluated from the 
condition of stationarity of the accretion process, adopting the rate of plasma penetration 
into the magnetic field of the neutron star equal to the rate of mass capture from the 
stellar wind, and, correspondingly, to the rate of mass accretion onto the stellar surface. 
Taking D e ff = £>b, where 

Db = «b/|^ (21) 
2eB(r m ) 



is the Bohm diffusion coefficient, and solving the equation M m (r mca ) = LxR n s/GM m , 
we find 

r mca ,2xl0 8 cmx a 2 ^ ^ ^ ^ (22) 

Here e is the electron charge, T(, — To/10 6 K is the plasma temperature and B(r m ) is the 
magnetic field strength in the magnetopause. Oo.i = «b/0.1 is the efficiency parameter 
normalized according to observational results and numerical simulations of solar wind 



penetration into the Earth magneto sphere( see y 81. 139] and references therein). 

Combining expressions (fT9b and (|22l) and solving inequality 2nl v S( j ^ K^(r mca ), we 
find 

|v sd | J> 2xl0- 14 Hzs- 1 ^ t% ? /13 m- 14 / 13 / 4 - 1 r 6 3/13 4g (23) 

-1 / d \ 17/13 
x ^57/13/ A ^ /AN 
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Thus, the observed spin-down rate of the neutron star in the X-ray pulsar 4U 2206+54 
can be explained in terms of magnetic accretion scenario assuming that the magnetic 
field strength on the neutron star surface is close to the value based on the analysis of 
X-ray spectra of this pulsar, i.e. ~ #crsf- 



4. CONCLUSIONS 

Our study of rapid spin-down of the neutron star in 4U 2206+54 together with presented 
earlier picture of spin evolution in another long -period X-ray pulsar GX 301-2 \ Tj^\ 



prove that application of the magnetic accretion scenario to the modeling of mass 
transfer in high-mass X-ray binaries is rather promising. One of the most important 
conditions for applicability of this scenario is magnetization of optical components in 
these systems which are usually giants and super-giants of early spectral type. Recent 
results of spectral and spectropolarimetric observations of massive stars (see, e.g. rf40l 
Hjl and references therein) show that strong magnetization is not unique among these 
objects. Magnetic fields strength in the range from hundreds to thousands Gauss (and in 
some cases in excess of 10 kG) have been measured in more than 20 early-type stars. The 
value of parameter /3 in the stellar wind from these objects is of the order of unity on the 
scale of orbital separation in close binary systems [42]. Another criterion for realization 
of magnetic accretion in a binary system is a moderate velocity of stellar wind satisfying 
the condition v cr < v re i < v mca . 

Thus, the systems with magnetic accretion occupy intermediate position between the 
systems with Keplerian accretion disks (v re i ^ v cr ) and the systems in which accretion 
onto a neutron star magnetosphere occurs in a form of free-falling quasi- spherical flow 
(v re i ^ v mca )- Finally, it should be noted that magnetic accretion scenario turns out to be 
the most effective in interpretation of long-period neutron stars whose magneto spheric 
radius is significantly less than the corotation radius. The spin-down mechanism based 



on matter outflow from the magneto spheric boundary (see, e.g. [|43|]) cannot be effective 
for these stars, while the spin-down torque caused by magnetic viscosity in the accretion 
flow remains significant. 
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